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Abstract 

It is shown that the K — > rnr matrix elements of the four-quark operator Q7, generated 
by the electroweak penguin-like diagrams of the Standard Model, can be calculated to 
first non-trivial order in the chiral expansion and in the 1/N C expansion. Although 
the resulting B factors B\ ' and B\ ' ' are found to depend only logarithmically on 
the matching scale /i, their actual numerical values turn out to be rather sensitive to 
the precise choice of /i in the GeV region. We compare our results to recent numerical 
evaluations from lattice-QCD and to other model estimates. 



1 Introduction 



In the Standard Model, the physics of non-leptonic -fT-decays is described by an effective 
Lagrangian which is the sum of four-quark operators Qi modulated by c-number coefficients 
q (Wilson coefficients) 



cl* 31 - 1 = -^v ud v: s J2ci(ii)(Km-\- h.c. . ( j ) 



where G F is the Fermi constant and V u d, V U s are the appropriate matrix elements of flavour 
mixing. This is the effective Lagrangian which results after integrating out the fields in 
the Standard Model with heavy masses (Z°,W ± ,t,b and c), in the presence of the strong 
interactions evaluated in perturbative QCD (pQCD) down to a scale \x below the mass of 
the charm quark M c . The scale /U has to be large enough for the pQCD evaluation of the 
coefficients c% to be valid and, therefore, it is much larger than the scale at which an effective 
Lagrangian description in terms of the Nambu-Goldstone degrees of freedom (K, tt and rj) of 
the spontaneous SU(3)l x SU(3)r symmetry breaking (SxSB) is appropriate. Furthermore, 
the evaluation of the coupling constants of the low-energy effective chiral Lagrangian which 
describes strangeness changing \AS\ = 1 transitions cannot be made within pQCD because at 
scales fi < 1 GeV we enter a regime where SxSB and confinement take place and the dynamics 
of QCD is then fully governed by non-perturbative phenomena. 

In this letter we shall be concerned with two of the four-quark operators in the Lagrangian 
in eq. @), the operators 

Q 7 = 6(Il7 m <Zl) ^RlnQR) ) ( 2 ) 

q=u,d,s 

and 

Q 8 = -12 ^2 e q(sLqn)(qRdL) , (3) 

q=u,d,s 

where e q denote quark charges in units of the electric charge, and summation over quark 
colour indices within brackets is understood. The operator Qi emerges at the Myy scale 
from considering the so-called electroweak penguin diagrams. In the presence of the strong 
interactions, the renormalization group evolution of Q7 from the scale Mw down to a scale 
fM < M c mixes, in particular, this operator with the four-quark density-density operator Qg. 
These two operators, times their corresponding Wilson coefficients, contribute to the lowest 
order O(p ) effective chiral Lagrangian which induces |A5*| = 1 transitions in the presence of 
electromagnetic interactions to order 0(a) and of virtual Z° exchange, i.e., the Lagrangian 

4T 1 = # kvo«) + *c . <4) 

Here, U is the matrix field which collects the octet of pseudoscalar Goldstone fields, Qr = 
diag[2/3, —1/3,-1/3] is the right-handed charge matrix associated with the electromagnetic 
couplings of the light quarks, and A^ 3 '* is the effective left-handed flavour matrix (A^ 23 ^ 
Si2$3j = 1,2,3). Under chiral rotations (Vl,Vr): 

U^V R UVl, Qr^V r Q r V ] r , X^^VlX^WI, (5) 

and the trace on the r.h.s. of eq. @ is an invariant. Actually, this is the only possible invariant 
which in the Standard Model can generate \AS\ = 1 transitions to orders O(a) and O(p ) in 
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the chiral expansion. With the given choice of the overall normalization factor in front of the 
r.h.s. of eq. (f|), the coupling constant h is dimensionless and, a priori, of order 0(iV c 2 ) in the 
1/N C expansion |2]-Q. This coupling constant plays a crucial role in the phenomenological 
analysis of radiative corrections to the K — > irir amplitudes ||. The determination of h is 
needed for a reliable estimate of these corrections Q The purpose of this note is to show that, 
following recent work reported in ref. fj\ , one can calculate the K — > tttt matrix elements of the 
Qi operator (and, therefore, the so called Bj 1 ' and Bj ' ' factors) to first non-trivial order 
in the chiral expansion and the 1/N C expansion This implies that (at least) the contribution 
to the constant h from the Qi and Qs terms in the four-quark effective Lagrangian can be 
calculated to first non-trivial order in the 1/N C expansion, a first step towards the required 
goal. The leading O(N^) contribution to h vanishes trivially. It could only come from the 
bosonization of the factorized Qs operator times its Wilson coefficient, but this coefficient is 
subleading at large N c 0. The next-to-leading contribution comes from the bosonization of 
the unfactorized Qj operator and its mixing via gluonic interactions with the Qs operator. It 
involves the same two-point function which governs the electroweak 7r + — 7r° mass difference 
in large--/V c QCD and this is why we are able to compute the contributions from the Qi 
and Qs terms of the four-quark Lagrangian to the constant h at the stated order in the 
1 /N c expansion. The bosonization of the factorized Qj operator and of the unfactorized Qs 
operator can only contribute to terms of order 0(p 2 ) (or higher) in the chiral expansion and 
they are, therefore, inoperational in the calculation of h. It turns out, however, that there is 
only a partial cancellation of the //-dependences generated by the product of the bosonization 
of the unfactorized Qj operator times its Wilson coefficient c-j with the //-dependence coming 
from the product of the Wilson coefficient eg, which is non-leading in the 1/N C expansion, 
times the bosonization of the factorized Qs operator. The full cancellation of //-dependences 
requires the consideration, as well, of the bosonization of other four-quark operators (in 
particular the unfactorized Q2 operator) in the presence of the electroweak interactions. This 
involves integrals of four-point functions which we have not yet fully explored within the 
framework of the 1/N C expansion. The constant h gets, therefore, contributions from other 
operators than just Q-j and Qs, and which we have not computed so far. We shall, therefore, 

(1/2) (3/2) 

concentrate here on the calculation of the factors Bj and Bj , and on their comparison 
to recent lattice QCD determinations [10-12], as well as to recent analytic determinations 
which have been made F] using the "effective action approach" of ref. [14]. 



2 Bosonization of Q7 and Qg 

As already mentioned, the bosonization of the operator Q7 is needed to next-to-leading order 
in the 1/N C expansion. The problem is entirely analogous to the bosonization of the operator 
Qlr = (qlI^QlQl) (qrI^QrQr) which governs the electroweak 7r + — ir° mass difference and 
which has been recently discussed in ref. @. Because of the LR structure, the factorized 
component of the operator Qj, which is leading in 1/N C , cannot contribute to order 0(p ) in 
the low-energy effective Lagrangian. The first non-trivial contribution from this operator is 

1 See ref. || for a recent discussion of the size of these corrections. 
2 A similar observation has also been made by J. Donoghuc 
3 See e.g. Buras's lectures Ej. 
4 See e.g. ref. H and references therein. 
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next -to-leading in the 1/N C expansion and is given by the integral [] 

A 

(27T) 



Qr - -sig,u J 70nf i? ( (? ) tr (ux^q r ) , (6) 



involving the two-point function 

n^(g) = 2i J d 4 x e**(0|T (^(x)^(O)t) |0) , (7) 

with currents 

^ = %(x) 7 ^(l- 7 5)?i(x) and R» = <?;(x) 7 ^(l + j 5 )qj(x) , (8) 

and z and j fixed flavour indices, i ^ j. This integral, which has to be evaluated in the chiral 
limit, where (Q 2 = — q 2 ) 

nf i? (g) = (gV-5 M V)n LR (g 2 ), (9) 

is divergent for large Q 2 and needs to be regulated. Before discussing this point, let us recall 
that in the large-iV c limit, the spectral function associated with Hlr,(Q 2 ) consists of the 
difference of an infinite number of narrow vector states and an infinite number of narrow 
axial-vector states, together with the Goldstone pion pole: 

-lmU LR (t) = ]T flM 2 b(t -M 2 )-Y, flM 2 A 5(t - M 2 A ) - f 2 5(t) . (10) 
71 V A 

Since Hlr(Q 2 ) obeys an unsubtracted dispersion relation, we find that 

Furthermore, in the chiral limit of QCD, the operator product expansion (OPE) applied to 
the correlation function Hlr(Q 2 ) implies 

lim Q 2 U LR (Q 2 )^0 and lim Q 4 U LR (Q 2 ) -> 0. (12) 

Q 2 — >oo Q 2 ^oo 



These relations result in the two Weinberg sum rules |ji 

E/v^v-E/i M i = /" and E/v<-E/i M l = °- ( 13 ) 

V A V A 

The usual prescription |16| for the evaluation of integrals such as (||) consists in taking a 
sharp cut-off in the (euclidian) integration over Q 2 , 

Qi - -6^2 £ d Q 2 Q 2 (-Q 2 n L dQ 2 )) tr (ux {23) u^q r ) . (u) 



5 Up to the replacement of Ql by A^ 23 ' , the bosonization of the operator Qi follows from the same procedure 
as described for the operator Qlr in ref. [Q. 
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With the constraints between the couplings and masses of the narrow states coming from the 
Weinberg sum rules (|l3|), the integral on the r.h.s. of eq. (|6j) becomes then only logarithmically 
dependent on the ultraviolet scale A, with the following result 



r A2 / 
J dQ 2 Q 2 {-Q 2 U LR (Q 2 



A 1V1 A V V 



(15) 



Notice that if only the contribution from the Goldstone pole had been taken into account, the 
resulting expression would have displayed a polynomial dependence on the cut-off scale A. 
Another possibility is to evaluate the integral (||) within a dimensional regularization scheme, 
say MS, in which case one obtains the same result (15), but with the correspondence between 



the cut-off A and the MS subtraction scale fi given by 

A = /i-e5. (16) 

On the other hand, at the level of approximation that we want to attain, the bosonization 
of the Qs operator is only required to leading order in the 1/N C expansion, because its Wilson 
coefficient is already subleading. To that order and to order 0(p ) in the chiral expansion 
it can be readily obtained from the bosonization of the factorized density currents, with the 
result [] 

Qs -> -12 Ub&J tr (uxftrfQn) , (17) 

where B is the low energy constant which describes the bilinear single flavour quark condensate 
in the chiral limit, B = —(ifiip)/f 2 . 



As discussed in ref. [18] there is a constraint that emerges in the large- N c limit which 
relates the leading d = 6 order parameter in the OPE of the Ulr two-point function [JH]] to 
couplings and masses of the narrow states: 

lim Q 6 U LR (Q 2 ) = -47r 2 (^ + 0(a 2 ))(^) 2 

= E fv M v - £ f\ M l ■ (18) 

V A 

This relation provides part of the cancellation between the \i dependence of the bosonization 
of the operator Qi{\i) with the short-distance dependence on \x in the Wilson coefficient of Qs', 
but, as already mentioned, and contrary to the simple case of the electroweak contribution to 
the 7r + — 7T° mass difference discussed in ref. 0], this cancellation here is incomplete. In this 
respect, we wish to comment on an important point concerning the scale dependence in the 
relation in eq. (|lg|). The term on the r.h.s. of the first line results from a lowest order pQCD 
calculation of the Wilson coefficient. Using the renormalization group improvement to one 
loop, this result becomes 

-^ m 2 ^-^f'(z^m\ h , (19) 

7T —Pi y 7T J 

with th e scale invariant quark condensate (the analog to invariant quark masses). In 

the large- N c limit, f3% — > ~^-N c and 71 — ► |iV c . The Q 2 dependence of the one-loop result 

6 See e.g. the lectures in ref. |17| and references therein. 



4 



/ 9 9 X -2 / 11 

is indeed rather mild ~ (logQ /A^g-J , but it does not go to a constant as the exact 
large-iV c result in the second line of eq. ( |lg| ) demands. This mismatch is due to the fact 
that pQCD is at best an approximation. It may happen that a two-loop renormalization 
group improvement of the OPE result approaches a constant behaviour at large Q 2 in a 
better way. In any case, this is a typical example of unavoidable mismatches that one will 
encounter between non-perturbative evaluations of matrix elements and short-distance pQCD 
evaluations, which are necessarily only approximate. In general, however, it is well known 
that at higher orders ambiguities will mix the short-distance coefficients of different powers 
in the OPE. It is difficult to imagine how to avoid these uncertainties in a final matching 
between short-distances and long-distances unless a breakthrough is made in understanding 
the relationship between pQCD and full QCD. 

Finally, if we restrict ourselves to the Lowest Meson Dominance (LMD) approximation 
to large- N c QCD discussed in ref. |2(| and identify My with M p , our calculation (in the MS 
scheme) gives the following contributions to the constant h coming from the operators Q7 
and Qs, 



^h[Q 7 ] = -180,0*) A 



log^-21og2 

p 



and ^h[Q 8 ]=-487r 2 c 8 ( M )^^, (20) 



where the relation (|l6|) is understood. 



3 The B Factors B^ 1/2) and B^ /2) 

The bosonic expression of Q7 given by eqs. ( |l4| ) and ( |l5| ) enables us to compute the K — > -kit 
matrix elements induced by this operator which, following the usual conventions, we express 
in terms of the following isospin amplitudes 

{Qi)i = {{^)i\Qi\K°), 1 = 0,2. (21) 

To leading order 0(jP) in the chiral expansion and to next-to-leading order in the 1/N C 
expansion, 0{\/ y/W^) for K — ► tttt amplitudes, we obtain the following result 



(Q 7 )o = V2(Q 7 ) 2 (>V ' :i 



16*2/2 



A 1V1 A V m V 



(22) 



It has become customary (rather unfortunately) to parameterize the results of weak matrix 
elements of four-quark operators Qi in terms of the factorized contributions from the so- 
called vacuum saturation approximation (VS A) , modulated by correction factors B^ 1 ^ , AI = 
1/2, 3/2. In the case of Q7 and Qs, one then has 

(Qr)^ = \X+ ±-{Z + AY) , (Q 7) vsa = -^X + , (23) 

and 

(Qs)^ = \{Z + m + , (Q 8 )^ SA = ^ Y ~^-X. (24) 
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The quantities X, Y and Z which appear in the above expressions are the same as those 
usually found in the literature, i.e., 



X 
Y 



sftU{M 2 K -Ml) + 0(p 4 



Ml 



m s + m di 



+ 0(p 2 ) , 



(25) 



There is no theoretical justification to consider the VSA as a good limit of any kind in QCD f\. 
This is reflected by the fact that e.g., the terms proportional to in eqs. and (pi]) do 
not correspond to the correct 1/N C expansion and, in fact, the contributions of Y and Z 
to (Q7)7 SA , although suppressed by a factor 1/N C , are actually numerically dominant, since 
(throughout, Mk denotes the neutral kaon mass) 



Y 



X 



Ml 



Ml 



Ml U s + m d 



11AX 



0.158 GeVV 



m s + m d 



■ 



(26) 



Furthermore, whereas X is scale independent, Y and Z depend on the MS subtraction scale 
\i through the quark mass term m s + m d in the denominator. We shall nevertheless follow 
these conventions, if only to be able to compare our results to those existing in the literature. 
The corresponding B factors are then defined as 



B 



(1/2) 



{Qi)o 
(Qi)F A 



B 



(3/2) 



(Q 



i/2 



(27) 



In the sequel, we quote our results in the MS scheme. 

Considering first the operator Q7, and restricting the sums in eq. (^) to the LMD ap- 
proximation discussed in ref. p3], the previous calculation leads to the results 



B, 



(1/2) 



(m) 



X 



X + jriZ + W) 



1 + 



My 



2vr 2 /2 Ml -Ml 



log 



A 2 

My 



2 log 2 



(28) 



and 



B, 



(3/2) 



(m) 



X 



1 



My 



A^flMl-Ml 



log 



A2_ 

My 



2 log 2 



If in these expressions we take the value My 
adopted in lattice calculations, 



M p , we obtain, at the scale fi 



(29) 
2 GeV usually 



B? /2) ( ) u = 2GeV) 



19.5 



1 + 23), Bf /2 \n = 2 GeV) 



■-a 

6.6 v 



11.5) 



(30) 



for, say, (m s + m d )(fi = 2 GeV) = 158 MeV, the conventional reference normalization used in 
eq. (^6|). These values are both positive and greater than unity. One should however notice 
that these numbers are rather sensitive to the choice of the scale \i and/or to the value assigned 
to My : for instance, at \i = 2My je~> ~ 1.3 GeV the contributions between square brackets 

7 Obviously, it would have been much more reasonable to normalize results to the large-JV c result in the 
lowest order of the chiral expansion. 
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on the r.h.s. of eqs. and ( p9|) vanish exactly. In order to illustrate these uncertainties, we 
show in Fig. 1 the variation of these B factors for a reasonable range of values of the scale \i. 
The area between the solid lines corresponds to the choice (m s +m ( i)(fj, = 2GeV) = 158 MeV, 
with A-rj-j varied between 300 MeV and 450 MeV. The area between the dashed lines reflects 
the same variation of Ajjg, but for the extreme low value {m s + m^)(^ = 2 GeV) = 100 MeV 
quoted in some lattice results [^]. Also shown is a magnification of the region corresponding 
to values of (i around 2 GeV, the reference scale at which the lattice results are usually given. 



We find that below fi < 1.3 GeV the Bj 1 ^ 2 ^ and Bj 0/ *' factors can even become negative, a 
result which disagrees, drastically, with the positive values quoted at a "matching scale" of 
0.8 GeV which are found in the constituent chiral quark model |l3|]. It is, however, not clear 
how this "matching scale" is related to the MS scale, fi, in QCD. 



(3/2) 





Fig. 1 The B factors, B^ 1/2) and Bf ,2) , 
lines correspond to (m s + md)(fJ> = 2 GeV) = 
value (m s + m d )(fJ- = 2 GeV) = 100 MeV. 



[L (GeV) 

as a function of the n scale in GeV. The solid 
- 158 MeV; the dashed lines to the extreme low 



Several matrix elements of four-quark operators have also been obtained in numerical 
simulations of lattice-QCD. These numerical evaluations, however, are based on a yet different 
definition of the B parameters. The lattice definition uses a current algebra relation between 
the K — » TTir and the K — > ir matrix elements which is in fact only valid at order O(p ) in the 
chiral expansion. Thus, in the case of Qj discussed here, the lattice-QCD B factor, which we 
shall distinguish with the symbol tilde on top, is defined by the ratio 



B 



(3/2) 



(* + \Q 



(3/2) 



(3/2) 



K- 



\VSA 

/o 



(31) 



where the matrix element in the denominator is evaluated in the chiral limit, as indicated by 
the subscript "0", and the operator Q-j has been decomposed into its AI =1/2 and AI = 3/2 
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components, Q-j = Q^ 2 ^ + Q 



,(3/2) 



with 



Q^, 3/2) = 2{s L ^d L )[u R 'y ll u R - d R <y^d R ] + 2(s L ^u L ){u R 'y^d R ) . 
The latest results obtained by various groups are consistent with each other: 



B!f /2) (fi = 2GeV) 



0.58 



+5+2 
-4-8 



(32) 



(33) 



latt | 0.61(11) 

The first value has been taken from ref. [10|, while the second one arises from the results of 



refs. [11] and |12j, translated into the MS scheme [[22]]. In the vacuum saturation approxima- 
tion which the lattice community uses, one obtains 



(<n + \Q? /2) \K^ 



VSA 



2 (Hp 
N c ft 



(34) 



where in the second expression, following common practice, we have traded the dependence 
with respect to the condensate for the dependence on the strange quark mass, using the Gell- 
Mann-Oakes-Renner relation, which holds in the chiral limit. The calculation based on the 
bosonized expression for this operator which we have discussed gives 



{7: + \Q ( f /2) \K + ) 



Svr 2 / 2 



V 



4vr 2 



My 21og2-log 



A 2 
M 2 



(35) 



where in the second line we have used the LMD approximation discussed in ref. |20|]. In this 
approximation and with the same numerical input as in (|30| ) , we obtain B^ 2 ^ 1 (/x = 2 GeV) ~ 
+1.5. However, as shown in Fig. 2, the value of this B factor is again very sensitive to the 
choice of the matching scale fi. The results obtained for the smaller value of the strange 
quark mass (i.e. for larger values of the condensate) are in agreement, within errors, with 
the numbers quoted in eq. ( p3| ) . Unfortunately, the values of m s (or of the condensate) and 
of the B factors obtained from the lattice are usually not quoted together. 

In the case of Q%, the large-iV c limit of eq. (p4|) simply reproduces the result obtained 
from the bosonized expression (|j~7|). Since subleading corrections in the 1/N C expansion of 



this operator are not yet available, further comparison with lattice results for B 
be postponed. 



(3/2) 



has to 
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1.8 1.9 2 2.1 2.2 



1 1.5 2 2.5 

(I(GeV) 

Fig. 2 T/ie -B factor, B^ 2 ^ [see eq. jSl\)], as a function of the scale in GeV. The solid 
lines correspond to (m s + m<i)(|U = 2 GeV) = 158 MeV; i/ie dashed lines to the extreme low 
value (m s + m d )(fj, = 2 GeV) = 100 MeV. 

4 Conclusions and Outlook 

The expressions in eqs. (|2~2j ) and (|35|) are the main results reported in this letter. They are a 
first step towards a systematic evaluation of weak matrix elements in the chiral expansion and 
to first non-trivial order in the 1/N C expansion. Our final aim, however, is to obtain values 
for the coupling constants of the low energy AS = 1 and AS = 2 chiral effective Lagrangian 
directly; i.e., constants like h in eq. (|j) and not of individual matrix elements of four-quark 
operators. It is encouraging from the results obtained so far, to find such simple analytic 
expressions which exhibit only a logarithmic dependence on the matching scale \x\ however, 
the fact that the numerical results are so sensitive to the choice of [i in the GeV region is 
perhaps an indication that one should be extremely cautious in the evaluation of errors of B 
factors, in general, both in model calculations and in lattice QCD numerical simulations. 
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